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ABSTRACT 


The focus of this research is to accurately model the tunnel junction interconnect 
within a multi-junction photovoltaic cell. A physically based 2-D model was created in 
Silvaco Inc.’s software to model the quantum tunneling effect that is realized within a 
multi-junction cell. The tunnel junction interconnect is a critical factor in the design of 
multi-junction photovoltaics and the successful modeling of the junction will lead to the 
ability to design more efficient solar cells. The quantum tunneling effect is based on the 
non-local band-to-band and trap assisted tunneling probability described by the Wentzel- 
Kramers-Brillouin (WKB) method. 
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EXECUTIVE SUMMARY 


The production of electrical power is of growing concern in today’s world. 
Renewable energy has been on the forefront of this issue and continues to provide 
promising results. Biofuels, wind turbines, hydroelectric, and geothermal energy are all 
sources of renewable energy. However, solar energy has proven to be a viable candidate 
for many niche applications. Whenever there is a need for a small, lightweight, low to no 
maintenance, completely independent source of energy, photovoltaic solar cells routinely 
provide a leading solution. 

Modern solar cell technology has been around since the 1950s and has greatly 
matured over the last half century. Solar power has developed into a multi-billion dollar 
a year industry [1] providing electric power to everything from spacecraft to automobiles. 
Although photovoltaic cells are a proven technology, there is always room for 
improvement. One major disadvantage of photovoltaic technology is the power density 
of the cells, which depends on the cell type and cell location. 

One solution to improve photovoltaic power density is to place cells of different 
material in tandem. These cells are connected with a tunnel junction in order to minimize 
shadowing of the sub-cells and maximize the number of photons absorbed. According to 
the DOE, a tunnel junction interconnect is a Tier 2 Technical Improvement Opportunity 
[1]. The potential for such a highly optimized photovoltaic cell depends on the ability to 
accurately model the quantum tunneling phenomenon in the tunnel junction. 

The focus of this research was to model the quantum tunneling phenomenon in a 

© 

degenerately doped p-n junction. This was achieved using Silvaco Inc.’s ATLAS 
software. Silicon (Si), gallium arsenide (GaAs) and indium gallium phosphide (InGaP) 
tunnel junctions were realized using this software and the affects of varying doping levels 
in response to a forward bias was investigated. 

The validity of the model was tested by implementing the successfully developed 
tunnel junction into a multi-junction solar cell with a known I-V curve. The multi¬ 
junction cell model was developed by Panayiotis Michalopoulos [2] and the model was 


XV 



compared to the experimental data found by the Central Research Laboratory, Japan 
Energy Corporation [3]. The simulated power density was calculated to be within 3.8% 
differenee from the experimental results. 
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I. 


INTRODUCTION 


A. BACKGROUND 

Photovoltaic solar cells are an exeellent source of renewable energy in countless 
applications. From Mars rovers to orbiting satellites, the need for a lightweight, 
renewable, low maintenance energy source is neeessary to extend equipment lifetimes 
and remain at the leading edge of technology and exploration. UAVs, sensor nodes and 
entire command posts are all using solar energy to power their systems to remain remote, 
untethered tools of our military. Even the modern war fighter is beeoming more and 
more reliant on gadget technology including night vision, satellite radios and GPS units, 
whieh are routinely issued to provide an advantage over an adversary. Considering these 
and many other photovoltaie applications, efficiency and size are paramount in researeh 
and development. 

Currently single eell photovoltaic energy sources are restricted to a small 
speetrum of emitted energy depending on the semiconductor used. Operating at between 
25 to 30% effieiency, the need for more power only inereases array weight by requiring 
more eells to produce a specified output power. To increase efficiency, multi-junction 
solar cells have been developed in order to absorb a wider array of the solar spectrum. 
By stacking cells of eompeting semiconductor material, efficiency can theoretically be 
inereased to upwards of 60%, whieh would lead to smaller arrays with increased power. 

The one flaw to simply staeking eells one on another is the parasitie junetion 
developed at the junetion between the solar cells. This problem is solved by 
manufacturing a tunnel junction between each cell in order to ensure that the eurrent flow 
in the eell and at the junction is in the same direetion. 

B. OBJECTIVES AND APPROACH 

The purpose of this thesis is to aecurately model the quantum tunneling 

phenomenon within a degenerately doped p-n junction. The junction will be modeled 

using Silvaco Data System Ine.’s ATLAS® software. ATLAS® is a multi-physics 

semiconductor simulation tool which provides the user with the ability to design and 
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evaluate ehanges to a produet before fabrieation. Onee designed, the developed junetion 
will be applied to a multi-junetion photovoltaie eell in order to inerease the effieieney 
while deereasing the development time and money spent. 

C. RELATED WORK 

This researeh is related to a Naval Postgraduate Sehool thesis eompleted in Mareh 
2002 by Panayiotis Miehalopoulos titled A Novel Approach for the Development and 
Optimization of State-of-the-art Photovoltaic Devices Using Silvaco. Although this 
thesis was very sueeessful in optimizing a multi-junetion photovoltaie eell, the forward 
bias I-V eharaeteristie of the tunnel junetion was not modeled. The ability to model this 
detail will further press the optimization of photovoltaie teehnology. 

D. ORGANIZATION 

The following two ehapters are a diseussion of semieonduetor physios and how 
they apply to a standard p-n junetion. Chapter IV extrapolates this information and 
applies it to a degenerately doped tunnel junetion. The simulation software along with 
the tunnel junetion simulation oode is disoussed in Chapter V. The results of the tunnel 
junetion model are presented in Chapter VI followed by the suooessfully implemented 
dual-junotion solar eell in Chapter VII. Finally, the reoommendations and oonolusions 
are presented in Chapter VIII. 
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II. SEMICONDUCTOR PHYSICS 


The basic material properties introduced in this section are essential to the 
understanding of the following research. Although there are more detailed works on this 
subject, this discussion will be limited to the essential properties related to silicon (Si), 
germanium (Ge), gallium arsenide (GaAs) and indium gallium phosphide (InGaP); 
semiconductors typical to photovoltaic applications. 

A. MATERIAL SCIENCE AND BOHR’S MODEL 

One common way to categorize a material is by its electrical properties. 
Depending on the level of an element’s resistivity, it can be categorized as an insulator, 
conductor, or semiconductor. The ability of a material to resist or conduct electricity is 
dependent on many factors, including lattice structure, free electrons, energy bandgap, 
and temperature. Some materials have very discrete electrical properties which define 
them as either an insulator or conductor. However, other materials such as silicon (Si) 
and gallium arsenide (GaAs) can act as either an insulator or conductor and are therefore 
considered semiconductors. Figure 1 shows the typical range of conductivities for 
insulators, conductors, and semiconductors. 
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Figure 1. Resistivity for various material types, from [4]. 


Bohr’s model can be derived as a first order approximation of quantum mechanics 
yet its simplicity is valuable to the understanding of the electrical characteristics of 
materials on a large scale. According to Niels Bohr, atoms are comprised of three 
subatomic particles: a negative charged electron, a positive charged proton, and a neutron 
with no charge. Compared to the electron, the high mass proton and neutron form what is 
known as the nucleus, which sits at the center of the model. The electron, which is about 
1,836 times smaller than the proton, revolves around the nucleus at a specific energy 
level known as an orbit or shell. The three subatomic particles are held together by the 
electrostatic force between the electron and proton similar to the way planets are held in 
orbit by gravity. Each shell must contain a set number of electrons before the outermost 
shell, known as the valence shell, can be occupied. 

Modeling the simplest of atomic structures, the hydrogen atom, Bohr assumes that 
an electron in orbit around the nucleus can take on only certain levels of angular 
momentum {nfi ). This quantization of the electron’s angular momentum is then related 

to the electron’s energy level by Equation 2.1 below, where m^is the mass of a free 
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electron, q is the electron charge, the permittivity of free space, n is the electron 


orbit and h is the reduced Plank’s constant 
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Eqn. 2.1 


As explained by Equation 2.1 and Eigure 2, it is clear that as the orbits increase, 
the energy at that particular orbit also increases by a factor of the orbit squared. It also 
follows that for an electron to move from an outer orbit to an inner orbit {E^-E ^), it 

must give off some of its energy, which it does in the form of light at a discrete 
wavelength (Eqn. 2.2). 



Eigure 2. Bohr’s model of allowable electron orbits in a hydrogen atom, from [5]. 
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Eqn. 2.2 
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As I will discuss later, this requirement of an eleetron to lose energy in order to 
move from an outer shell to an inner shell is the same principle a solar cell uses to absorb 
light and generate power. 

B. CRYSTAL STRUCTURE 

The spatial arrangement of the atoms within a physical structure greatly affects 
the precise properties of a given material. As Figure 3 illustrates, a solid ean be elassilfed 
into three separate degrees of atomie order: amorphous, polyerystalline and crystalline. 
Amorphous structures are arranged in such a way that any one section within the 
structure is different from any other section. Crystalline structures are arranged such that 
the atoms are arranged in the same manner throughout the solid. Polyerystalline 
structures area a eombination of both amorphous and erystalline. Although it is possible 
for a semiconductor to be arranged in any of the three atomic order classifications, 
polyerystalline solids are most prominent in today’s market due to cost, ease of 
manufaeturing and adequate reliability in photovoltaie applications. 
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Figure 3. Degree of atomic order classification within solids: (a) amorphous (b) 

polyerystalline (c) crystalline, from [5]. 
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The atomic structure of a semiconductor can be further described by its lattice 
structure within a unit cell. The typical way to describe an atomic structure is to use the 
Brevais lattice system, which describes the symmetry of a crystal by the three- 
dimensional (3-D) geometric arrangement of the atoms in it. Although there are seven 
separate categories within the Brevais system with fourteen total arrangements, the 
simplest of the arrangements can be divided into three categories: the simple cubic, body- 
centered cubic and face-centered cubic (Figure 4). 



Figure 4. Brevais lattice cubics, from [6]. 


Taking a closer look at the lattice structure of an element with four valence 
electrons, such as carbon, silicon and germanium, there are small differences from the 
lattice structures described above and those of the group IV elements (those with four 
valence electrons). Figure 5 is a diamond lattice in which two face centered cubic lattices 
are interpenetrating. The blue and green spheres represent the first interpenetrating face 
centered cubic and the light blue sphere represents the second interpenetrating face 
centered cubic. The second face centered cubic is diagonally shifted Vi of the unit 
distance. All element IV semiconductors are described as having cubic-diamond lattice 
and are so named because each falls under the carbon element, which is the basic 
building block of diamond. Within this lattice structure each atom is equidistant from its 
four nearest neighbors. Silicon’s and germanium’s lattice constants are 5.437 A and 
5.660 A respectively. 
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Figure 5. Cubic-Diamond lattice structure, from [6]. 


C. CARRIERS 

Carriers within a material provide an avenue for the transport of an electric charge 
and therefore electric current. Similar to the Bohr model of the hydrogen atom, the 
electrons within a Si atom only occur at certain energy levels. However, since the silicon 
atom is a multi-electron atom, the associated shell model is slightly more complex. The 
innermost filled shells along with the nucleus are often referred to as the core of the atom 
because they are tightly bound by the electrostatic force between these shells and the 
nucleus. The remaining electrons of the atom that fill the outer shell have a weak bond to 
the core, and are referred to as the valence electrons due to their strong participation in 
chemical reactions. 


8 



















Silicon has fourteen eleetrons, ten of whieh oeeupy the two inner most shells 
(atom eore) leaving four valence eleetrons. When the lattiee is formed between the 
different atoms, there is a finite probability that at any temperature above zero Kelvin 
(OK), an eleetron from the valenee band will be knoeked free and an eleetron deficieney 
known as a hole will be left behind. The hole left behind retains a positive eharge equal 
to the absolute value of the eleetron eharge. The free eleetron along with the hole is 
known as an eleetron hole pair (EHP) and ean be deseribed by the phenomenon of 
ionization or generation. When a free eleetron and hole eombine they disappear, a 
proeess whieh is known as reeombination. The rate of reeombination is proportional to 
the number of free eleetrons and existing holes. When in thermal equilibrium the 
generation and reeombination keep the number of EHPs eonstant. Germanium reaets in 
mueh the same way although it has thirty-two total eleetrons, twenty-eight of whieh 
eonstitute the atom eore, leaving four valenee eleetrons. Reeombination and generation 
will be diseussed in more detail later. 

An intrinsie material has no impurities (doping) in it and will not produee EHPs at 
OK. However, the EHPs will inerease logarithmieally as a function of temperature 
(Eigure 6) and therefore eonduetivity will also inerease. 
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Temperature [°K] 


Figure 6. Intrinsic carrier density vs. Temperature, from [7]. 


D. ENERGY BANDS 

The energy band model is a useful way to visualize the differences between an 
insulator, conductor and semiconductor (Figure 7). When looking at a material rather 
than a free atom, the electrons within the material form energy bands known as a valence 
band and conduction band at a specific energy state. The difference between the valence 
and conduction band is known as the bandgap. The conductivity of a material is 
dependent on the number of free electrons in the conduction band. 

An insulator has a large bandgap that an electron must overcome in order for it to 
flow in the conduction band. Conversely, a conductor is described as material with the 
valence and conduction bands overlapping in which case the free electrons do not require 
any energy to traverse an energy barrier. As expected, a semiconductor has a bandgap 
that is not overlapping yet the energy difference between the valence and conduction 
band is minimal. 
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Figure 7. Energy band diagrams, from [8] 


Within the energy band there is a distinetion between a direet bandgap and an 
indireet bandgap based on the relative momentum between the valance band hole and the 
conduction band electron. A direct bandgap semiconductor is one in which the maximum 
energy level of the valence band and the minimum energy of the conduction band occur 
at the same momentum. Examples of direct bandgap semiconductors are GaAs and InP 
which are sometimes referred to as radiant recombination semiconductors because of the 
band-to-band recombination and release of energy in the form of photons or radiation 
wavelengths of light. 

Indirect bandgap semiconductors have maximum and minimum energy band 
levels at different momentums. Examples of indirect bandgap semiconductors are Si, Ge 
and GaP. The energy released during indirect recombination cannot be direct band-to- 
band and requires the assistance of a third party or defect in the material. This energy is 
given off in the form of phonons or quasi-particles of quantized sound waves. The 
energy band structure for Ge, Si and GaAs are compared in Eigure 8. This relationship 
will become more clear when discussing recombination and generation. 
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WAVE VECTOR 


Figure 8. Energy band strueture vs. wave veetor for Ge, Si and GaAs, from [7]. 


E. INTRINSIC SEMICONDUCTORS AND DOPING 

An intrinsie semieonduetor is one with no impurities in it. Therefore, it exhibits 
properties inherent to the pure semieonduetor sample whieh ean be used to identify the 
material. Under equilibrium eonditions, an intrinsie semieonduetor has an equal number 
of holes (p) and eleetrons (n) whieh is ultimately equal to the intrinsie earrier 
eoneentration (n^). The number of p and n are equal in an intrinsie semieonduetor 
beeause earriers are only ereated in pairs in a highly pure material. If an intrinsie 
semieonduetor bond is broken, a free eleetron and a hole are ereated simultaneously. 
Similarly, if an eleetron is exeited from the valenee band to the eonduetion band then a 
hole will simultaneously be ereated in the valenee band. A list of eommon 
semieonduetor materials and their intrinsie earrier eoneentrations at room temperature are 
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given in Table 1. Although these eoneentrations may seem high, they are only a small 
fraetion of the total number of bonds in the respeetive materials. 


Table 1. Intrinsie earner eoneentration at room temperature (300K). 


Carrier Concentration ( n .) per cm^ 

Material 

2x\0^ 

GaAs 

1x10'“ 

Si 

2xl0‘“ 

Ge 


Adding impurity atoms to a material changes its electrical properties, a process 
known as doping. Doping is used to control the number of electrons or holes in a 
material by the addition of a controlled amount of a specific impurity atom. The doping 
of intrinsic material results in either a p-type or an n-type material, meaning it has an 
excess number of holes or electrons respectively. With respect to IV elements, acceptors 
(III elements) are added to produce a p-type material and donors (V elements) are added 
to produce n-type materials. This is due to the number of valence electrons in the dopant 
atom. From Section II.A we know that IV elements have four valence electrons. If an 
element with three valence electrons is added to the diamond lattice the three valence 
electrons will not be able to complete all of the semiconductor bonds with the element IV 
atoms. Therefore, the excess electron will continue to be accepted by the dopant atom, 
and an excess hole will be allowed to wander throughout the structure without the 
creation of a free electron. Conversely, if a donor is added to the intrinsic element IV 
material then all the semiconductor bonds will be made, leaving an excess electron to 
roam throughout the lattice. Only a small amount of dopant, on the order of one atom for 
every 10^'’ atoms, is required to significantly affect the electrical properties of a material. 


13 




F. FERMI LEVEL 

The Fermi level describes the top of the collection of electron energy levels at OK 
and is based on the Fermi-Dirac probability equation (Eqn. 2.3). The Fermi-Dirac 
probability equation is a particle statistic function that determines the distribution of 
fermions (particles with half-integer spin) at a particular energy state at absolute zero. 


f(E) = 


1 

^(E-E,)lkT 


Eqn. 2.3 


At OK all the available energy states below will be filled by fermions of one 

and only one particle which are constrained by the Pauli Exclusion Principle. At low 
temperatures, any energy below the Eermi level will have a probability of 1 and any 
energy above the Eermi level will have a probability of essentially 0. As the temperature 
rises, EHPs are created and electrons of higher energy than the Eermi energy emerge in 
the conduction band. Eigure 9 illustrates this distribution for increased temperatures in an 
intrinsic material. 
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Eigure 9. Eermi distribution function, from [5]. 


The Eermi level can be manipulated by doping the material either with an 
acceptor or donor, resulting in a p-type or n-type semiconductor respectively. By raising 
or lowering the Eermi level, one can specify the required energy for an electron to 
transition from the valance band to the conduction band. As in Eigure 10, an intrinsic 
material has a Eermi level at the center of the band gap. An n-type material shifts the 
Eermi level closer to the conduction band producing shallow donors. Similarly, a p-type 
material shifts the Eermi level closer to the valence band producing shallow acceptors. 
With the Eermi level closer to the valence and conduction bands, only a small amount of 
energy is required to ionize a hole or electron into their respective bands. 
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Figure 10. Schematic band diagram, density of states, Fermi-Dirac distribution and the 
carrier concentrations for intrinsic, n-type and p-type semiconductors at thermal 

equilibrium, from [7]. 
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When the material is not under equilibrium eonditions the Fermi level is 
influeneed ereating Quasi-Fermi energy levels {Ep^ andF’^^) deseribed by the 

distributions in Equation 2.4 and Equation 2.5. It is easy to see from these two equations 
that at equilibrium Ep = Ep^ = Ep^ . 


fn(E) = 


1 

^{E-E,„)/kT 


Eqn. 2.4 


f,iE) = 


1 


^(E-E,„)lkT 


+ 1 


Eqn. 2.5 


As I will diseuss later with respeet to the tunneling effeet, it is also possible to 
degenerately dope the material in order to push the Eermi level into the valenee or 
eonduetion band. However, it is important to note that a heavily doped material may 
exhibit unintended eleetrieal properties depending on the distribution of the dopant. 
Beeause the distribution of the dopant may not be homogenous throughout, different 
energy fluetuations may be found at different positions within a strongly eompensated 
semieonductor. This is due to the mismateh of the lattiees between the two materials, 
whieh ean greatly affeet the eonduetivity of a material. 

G. MOBILITY 

Mobility (p) is the ratio of average speed (drift speed) to the applied field of a 
material and is analogous to permittivity ). Ideally (in a vaeuum) an eleetron inside of 

an eleetrie field will eontinuously aeeelerate. However, in a material the eleetron will 
aeeelerate only until it eollides with an atom in the lattiee. As expeeted, the higher 
doping in a material will result in a higher probability of a eollision with atoms in the 
lattiee and therefore mobility will deerease. Eigure 11 illustrates this point for some 
eommon semieonduetors. 
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Figure 11. Mobility vs. doping level, from [7]. 


H. CURRENT DENSITY 

The total particle current in a material is equal to the total current density (J) of 
both the holes (Jp) and electrons () which results from the drift and diffusion current 
densities. Diffusion current (I^,) is the result of electrons or holes in a concentrated area 
spreading themselves evenly though a material (Figure 12) in order to achieve a zero 
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concentration gradient (Vn = 0 for holes, Vp = 0 for eleetrons) throughout. Diffusion 
eurrent density (for holes, Jj^|jjj.pfor electrons) is analogous to the majority earrier 

concentration (doping) of a material and is deseribed by Equation 2.6 and Equation 2.7 
where Dp and are the eonstants of proportionality. 


Id 

-► 



Jp|d.ff = -qDj,yp Eqn. 2.6 

JNidiff = Eqn. 2.7 


Drift eurrent (Ig) is the result of the introduction of an electric or magnetic field 

and is analogous to the minority carrier eoncentration and is proportional to the intensity 
of the field (rf). Eigure 13 illustrates drift eurrent in the presence of an eleetric field and 
Equation 2.8 and Equation 2.9 describe the mathematical model of drift eurrent as a 
function of mobility (p) and applied electric field. 
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Figure 13. Drift current in the presence of an electric field, from [2]. 

Eqn. 2.8 
Eqn. 2.9 


Jpldrift 
■^Nldrift ^ 


Is 

► 



Both the electron (//^) and hole (//^) current densities depend on their associated 

diffusion and drift currents (Eqn. 2.10 and Eqn. 2.11) which ultimately lead to the total 
current density (Eqn. 2.12). 


JP = Jp|drift + Jp|diff = 9/^pE^ - 

Eqn. 2.10 

■^N “ drift -^N diff “ qPrP^' ~ q^^^ 

Eqn. 2.11 

J — + Jp 

Eqn. 2.12 


1. GENERATION AND RECOMBINATION 


Generation is the process whereby electrons and holes (carriers) are created. 
Conversely, recombination is the process whereby electrons and holes are annihilated. 
Both generation and recombination are order-restoring mechanisms to return the system 
to its equilibrium state (p=n= n .) after the carrier excess or deficit in the semiconductor is 
stabilized or eliminated. 
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There are three important recombination-generation methods as shown in Figure 
14. Band-to-band recombination (a) is the direct annihilation of a valance band hole and 
a conduction band electron when the hole and electron stray into the same vicinity within 
the semiconductor lattice. The loss of energy experienced by the electron must release an 
energy quantum equal to the total energy lost. This loss is typically realized in the form 
of light as described in Equation 2.2. 

The R-G (recombination-generation) center recombination is caused by defects or 
special impurity atoms within the lattice and takes place only in special locations within 
the semiconductor known as R-G centers. The R-G center concentration is typically very 
low compared to the donor and acceptor concentration yet leads to the introduction of 
various energy levels near the center of the band gap known as the trap energy {Ej.). 

The R-G center recombination is a two-step process and can take place in one of two 
ways as represented in Figure 14(b). The first method is when a carrier (hole or electron) 
is trapped at the impurity location until the opposite carrier is attracted to the first and 
recombination annihilates them both. The second is when an electron falls into the 
impurity location where it will initially lose energy equal to Once a hole 

comes into vicinity, the electron is attracted into the valence band after which it will lose 
the rest of its energy equal to \Ej - Ey | before it recombines with the hole and disappears. 

This type of recombination will cause quantum lattice vibrations which will release 
thermal energy and ultimately increase the temperature of the semiconductor. 

The third process is called Auger recombination (c), which is also a band-to-band 
recombination. However, this process occurs simultaneously with a collision between 
two like carriers. The energy transferred between the colliding carriers forces one carrier 
into the alternate energy band while the second carrier remains in the original band. The 
second carrier becomes highly energetic and then returns to its original energy level in 
small heat producing steps as it collides with the semiconductor lattice. This unique 
stepwise energy loss or “staircase” is defined as thermalization and also leads to an 
increased semiconductor temperature. 
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Figure 14. Energy band visualization of recombination and generation processes, 

from [5]. 


The generation process is simply the reverse of the recombination process. For 
band-to-band generation (d), thermal energy or light energy >£^g) absorbed by a 

semiconductor will excite electrons from the valence band into the conduction band 
creating an EHP. Similarly, for R-G center generation (e), thermal energy absorbed by 
the semiconductor will excite the electron enough to achieve the trap energy at the 
impurity location and then again to the conduction band. Carrier generation via impact 
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ionization (f) is analogous to Auger recombination, because the creation of EHPs is the 
result of energy released when a highly energetic carrier collides with the semiconductor 
lattice. 

J. SUMMARY 

The basic understanding of semiconductor physics leads to a more advanced 
discussion of a p-n junction, which is discussed in the next chapter. 
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III. P-N JUNCTIONS 


A. INTRODUCTION TO JUNCTIONS 

The successful implementation of a P-N junction is critical in modern electronics 
and semiconductor devices, including bipolar junction and metal on oxide transistors 
(BIT and MOSFET), rectifier and Zener diodes, and photovoltaics to name a few. 
Through the doping process, a P-N junction is formed to effectively control the electrical 
characteristics of a material. 


B. BASIC STRUCTURE 

A p-n junction is formed by a single crystal semiconductor where one side is 
doped with acceptors (p) and the other side is doped with donors (n). At the location of 
contact the metallurgical junction is formed as shown in Figure 15. Initially there is a 
very large majority carrier density gradient at the metallurgical junction. Equilibrium 
will be sought as the majority carriers in the p-type material (holes) will begin to diffuse 
into the n-type material while the majority carriers from the n-type material (electrons) 
will begin to diffuse into the p-type material. 



\ 


Metallurgical 

junction 




Hole diffusion 





Electron diffusion 

x = 0 


(a) (b) 

Figure 15. (a) Simplified P-N junction; (b) doping profde, from [9]. 
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The diffusing carriers from both regions will combine with the minority carriers 
on the opposite side and disappear. Assuming that there is no outside electric field 
connected to the junction, this diffusion can not continue indefinitely. The holes 
diffusing from the p-region will leave behind negatively charged acceptor atoms and the 
electrons diffusing from the n-region will leave behind positively charged donor atoms. 
The net charges in the p- and n-regions will induce an electric field in the metallurgical 
junction, which will form what is known as a space charge region or depletion region as 
shown in Figure 16. 
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Figure 16. P-N junction depletion region formed by diffusion, from [9]. 


When in thermal equilibrium, the Fermi levels on both sides of the junction are 
equal (Figure 17). The potential barrier known as a built in voltage (I 4 ,.) is responsible 
for maintaining equilibrium between the majority carrier electrons in the n-region and the 
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minority carrier electrons in the p-region. Similarly for the holes in the valence band, the 
built in voltage maintains equilibrium between the majority carrier holes in the p-region 
and the minority carrier holes in the n-region. As expected, the built in voltage depends 
on doping levels and A^) within the n- and p-regions (Eqn. 3.1). 
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Figure 17. P-N junction energy band diagram with Fermi level. 


Ki =—In 

q 


'Ma" 


Eqn. 3.1 


The width of the depletion region (W) is also related to the relative doping level 
of the two regions. This characteristic will become more important with the 
implementation of the tunnel junction (Section IV). The width of the depletion region (in 
cm) can be calculated by using Equation 3.2. 


1T = 




i ^ 

_ Jj 


Eqn. 3.2 


27 
























Applying different bias conditions between the two regions will disturb the p-n 
junction from equilibrium and the Fermi levels will no longer remain constant. For a 
standard doped junction, the three conditions of interest are reverse bias, forward bias and 
the breakdown voltage. 


C. REVERSE BIAS 

Applying a negative voltage (I 4 ) to the p-region with respect to the n-region 

results in a reverse bias condition as shown in Figure 18. The voltage applied attracts the 
majority carriers from their respective regions and leaves minority carriers in its place. 
This increases the voltage barrier and widens the depletion region. The difference 
between the two Fermi levels is equal to the bias voltage. Therefore, the total barrier 
voltage ) realized by the system is the sum of the built in voltage and bias voltage 
(Eqn. 3.3). This barrier is easiest to visualize by the use of an energy band diagram in 
Figure 19. 



® Hole 49 Elechou lecouibiiiiiig with a hole I I Neiihal charge 

9 Elech ou ^9 Hole recoiiibiiiiug with an elech ou | | Posih^•e charge 

. . Elechostatic field I I Negative charge 


Figure 18. Reversed bias p-n junction schematic, from [2]. 
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Eqn. 3.3 
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Eigure 19. Reversed bias p-n junetion energy band diagram. 


The depletion region width is given by Equation 3.4. When in reverse bias the 
probability that an electron will gain enough energy to overcome the potential barrier is 
very small and the junction will act like an insulator, not allowing current to flow. With 
the increased depletion region, diffusion current {Id) is very small and the total current in 
the system is approximately equal to the drift current {Is) (Eqn. 3.5). Ideally, this current 
would be equal to zero, but, in reality, there is a small amount of drift current in the 
system, as illustrated in Eigure 20. 


W = 


\2e,(V„+V,) 


[ ^ 

L NaN, Jj 


Eqn. 3.4 
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ho,ai =h+^D~h (very small) 


Eqn. 3.5 



Figure 20. Reversed bias p-n junction I-V characteristic. 

D. FORWARD BIAS 

Applying a positive voltage (1/) to the p-region with respect to the n-region is 

considered forward bias (Figure 21). The forward bias opposes the built in field 
developed in the depletion region and the region becomes narrower. When the junction 
is forward bias, the Fermi level in the p-region is increased and the Fermi level in the n- 
region is decreased to a point where electrons can overcome the small potential barrier 
(Eqn. 3.6), as illustrated by Figure 22. 
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Figure 21. 


Forward bias p-n junction schematic, from [2], 


(V„-VJ Eqn.3.6 

hotai=h+lD~lD Eqn. 3.7 
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Figure 22. Forward bias p-n junction energy band diagram. 


As the continuously supplied majority carriers from the n-region transition into 
the p-region, they become minority carriers where they diffuse and eventually combine 
with the majority holes. Similarly, the continuously generated holes in the p-region will 
be available to combine with the electrons traversing the potential barrier. Just as before, 
some of the current being produced is due to drift (Is) and some is due to diffusion (Id). 
However, in the forward bias condition the diffusion current is much greater than the drift 
current (Eqn. 3.7) due to the decreased depletion region. Figure 23 shows the I-V 
characteristic of a p-n junction in the forward bias condition. Notice that the current does 
not begin to flow until sufficient voltage is applied to overcome the potential barrier. 
This happens at approximately 0.7V for Si and O.IV for Ge. 
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Figure 23. Forward bias p-n junction I-V characteristic. 


E. BREAKDOWN/ZENER TUNNELING 

Adding reverse bias beyond a certain threshold voltage will result in the 
breakdown of the p-n junction. This breakdown phenomenon is due to either a Zener 
effect or an avalanche effect. 

The Zener or tunneling effect is similar to the tunneling which I will discuss at 
length in the next section. The difference here is that the Zener effect takes place when a 
minority carrier in the p-region valence band is able to penetrate through the energy 
barrier to an empty state in the n-region conduction band unaffected. This can only occur 
once the valence band of the p-region increases beyond the conduction band of the n- 
region due to the reverse bias. The exchange of minority carriers creates a large change 
in reverse current with a small change in reverse bias. Zener diodes are designed to 
achieve this result (Figure 24). 
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Figure 24. Zener breakdown p-n junetion I-V eharacteristic. 


The avalanche effect is also a minority carrier effect in which an electron from the 
p-region with a very large kinetic energy collides with atoms and breaks the covalent 
bonds. Often the EHPs created collide with other atoms to create even more EHPs, 
creating an avalanche-type effect. As more and more minority carriers are created, they 
are swept to the opposite region which creates a large reverse current. This runaway 
effect often leads to excess heat from the lattice vibrations and results in a damaged p-n 
junction. 

F. LOSSES 

Ideally there would be no losses in the system and the slope of the I-V curve 
would be perpendicular to the x-axis when on. However, since the semiconductor 
material has an inherent ohmic resistance due to atomic masses and impurities, there is a 
slight ohmic loss in the junction. Putting three regions together, a typical p-n junction I- 
V curve is realized (Eigure 25). 
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Figure 25. p-n junction I-V characteristic. 


G. SUMMARY 

Controlling the doping levels in the respective regions will shift the Fermi levels 
and affect the ability of an electron or hole to move through the lattice. Increasing the 
doping levels to a point where the Fermi levels are within the energy bands of the n- and 
p- regions is essential to tunneling phenomenon and will be discussed in the following 
Chapter. 
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IV. TUNNEL JUNCTION 


A. TUNNELING THEORY 

The tunneling phenomenon is a majority earner effect in which both the n- and p- 
regions are degenerately doped (very heavily doped with impurities) allowing an electron 
to tunnel through a potential barrier that, under normal conditions, would not be possible 
to penetrate. The tunneling is based on a quantum transition probability (7^) per unit 
time and is proportional to the exponent of the average value of momentum (^ (o) ) within 
the tunneling path across the depletion region (Eqn. 4.1). This momentum corresponds to 
an incident carrier with a transverse momentum of zero and energy equal to the Fermi 
energy. The tunneling probability equation will be further explained later in this section. 


= exp 


-2| " A:(x) dx 


Eqn. 4.1 


B. PURPOSE 

To achieve a higher power density, multiple cells with different I-V curves and 
decreasing bandgaps are placed in tandem as illustrated in Figure 26. As a result the 
incident light will be absorbed by the top cell allowing the unabsorbed photons to pass 
through to the sub cells. This results in the sum of the voltage produced by the individual 
cells. 


37 




Figure 26. 


Tandem solar cells. 


Unfortunately, stacking the cells in tandem would result in a parasitic junction 
formed by the two dissimilar regions of the different cells (Figure 27) and must be 
separated in order to allow current flow. 
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Figure 27. Tandem solar cell with parasitic junction. 
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One option is to mechanically stack the cells using an oxide/conductor contact as 
shown in Figure 28. However, the added contact material shadows the sub cells from the 
unused photons being passed through the top cell, limiting any efficiency gained by 
placing the cells in tandem. 
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Figure 28. Tandem solar cell example with oxide/conductor contact separation, from [2]. 

Another solution to this issue is to use a tunnel junction to separate the individual 
cells in the tandem cell. The tunnel junction will allow photons to pass through as well as 
pass current with minimal voltage loss (Figure 29). 
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Figure 29. Tandem cell with tunnel junction. 


C. THE JUNCTION 

A tunnel junction is simply a tunnel diode or degenerately doped p-n junction 
commonly used in switching and oscillator circuit applications. The difference between 
the tu nn el junction and the tunnel diode is that the tu nn el junction is grown as a single 
crystal joining two solar cells, whereas a tunnel diode operates as a discrete component 
within a circuit. The purpose of mentioning the tu nn el diode is for completeness, as 
tunnel diodes have been studied and used in electronics since their invention by Leo 
Esaki in 1958, who received the Nobel Prize in Physics in 1973 for his work. From this 
point on, there will be no distinction made between the two and I will focus on the tunnel 
junction and its application to photovoltaics. 

A tunnel junction is formed with degenerately doped regions in order to push the 
respective Fermi levels into the valence and conduction bands (Figure 30). As a result, 
the p-region valence band will contain holes and the n-region conduction band will 
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contain electrons at thermal equilibrium. The depletion region width is also affeeted by 
the degenerately doped regions compressing it to <100A, whieh is essential to the 
tunneling effeet and is signifieantly smaller than a eonventional p-n junetion. 


E 

Ec 

Ev 

Ef 


luglily p-doped 2 luglily n-doped 
semiconductor ^ semiconductor 


Figure 30. Degenerately doped tunnel junction energy band diagram. 


In order for quantum tunneling to occur, four criteria must be met. The neeessary 
conditions for tunneling are: (1) oecupied energy states exist on the side from whieh the 
electron will tunnel; (2) unoccupied energy states exist at the same energy levels as those 
electrons in the occupied energy states; (3) the tunneling barrier width is small and the 
potential barrier height is low enough that there is a finite tunneling probability; and (4) 
the momentum is eonserved in the tunneling proeess [2]. 

D. THERMAL EQUILIBRIUM 

As in the eonventionally doped p-n junetion, the two Fermi levels of the p- and n- 
regions are equal at thermal equilibrium. Figure 31 illustrates the tunneling junetion 

energy bands at thermal equilibrium as well as the eorresponding I-V plot. The red dot 
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on the I-V plot corresponds to the location of interest on the curve. Since both regions 
are doped to the point where the Fermi levels exist inside the valence and conduction 
band and are therefore considered degenerately doped. 



Figure 31. Tunnel junction energy band and I-V characteristic at thermal equilibrium, 

from [7]. 


E. REVERSE BIAS 

Similar to the Zener effect discussed in the last section, in reverse bias the 
minority carrier in the p-region valence band is able to tunnel through the potential 
barrier to the conduction band of the n-region. Unlike the Zener effect in a standard p-n 
junction, the reverse bias tunneling in this junction occurs as soon as a reverse bias is 
applied to the junction. The reason for the immediate response is the degenerate doping 
of the regions. Since the Fermi levels are inside the respective valence and conduction 
bands at equilibrium, there are empty states available for the minority carrier to fill as 
soon as reverse bias is applied. Therefore, a large reverse current is realized for a small 
change in voltage. Figure 32 illustrates the minority carrier electron tunneling to the 
conduction band of the n-region as well as the location on the I-V curve where reverse 
bias tunneling takes place. 
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Figure 32. Tunnel junction energy band and I-V characteristic in forward bias 

(reverse bias), from [7]. 


F. FORWARD BIAS 

At forward bias (Figure 33), the tunnel junction has available energy states in the 
p-region corresponding to filled energy states in the n-region. As the voltage is increased 
from zero to , the valence band energy of the p-region and the conduction band energy 

of the n-region become closer in value. 
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Figure 33. Tunnel junction energy band and I-V characteristic in forward bias 

(peak tunneling), from [7]. 


G. NEGATIVE RESISTANCE 

From Vp to V^,, tunneling is still occurring but fewer and fewer empty states are 
available. With a small increase in forward bias the number of electrons in the n- 

region directly opposite the empty states in the p-region continues to decrease which will 
result in a decreased tunneling current. The region from 1/^ to V,, is described as 

negative resistance (Figure 34) and will continue until the valence band of the p-region is 
equal to the conduction band of the n-region at which time the tunneling current will be 
zero. 
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Figure 34. Tunnel junction energy band and I-V characteristic in forward bias 

(negative resistance), from [7]. 


H. THERMAL CURRENT FLOW 

Increasing the forward bias even further beyond V,, will eventually result in the 

normal current flow of a p-n junction. As illustrated in Figure 35, the energy band 
diagrams of the p- and n-regions are reminiscent of the standard p-n junction with an 
applied forward bias. The relationship between a standard p-n junction and a 
degenerately doped tunnel junction can be further understood by comparing the two 
curves in Figure 36. As additional forward bias is applied, the current flow of the 
standard p-n junction and that of the tunnel junction merge. 
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Figure 35. Tunnel junction energy band and I-V characteristic in forward bias (thermal 

current flow), from [7]. 



Figure 36. Characteristic I-V curve of a tunnel (red) and conventional (blue) p-n junction. 
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I. DIRECT AND INDIRECT TUNNELING 

Just as there is a direct and indirect materials relationship, so too is there a direct 
and indirect tunneling process. Superimposing the E-k relationship at the classical 
turning points onto the E-x relationship of the tunnel junction reveals the direct and 
indirect tunneling differences. In Eigure 37 (a) and (b) the electrons in the vicinity of the 
minimum of the conduction band energy-momentum can tunnel to a corresponding value 
of momentum in the vicinity of the valence band. However, the distinction lies in the 
electron and hole momentum. 

Eor direct band tunneling to occur, the momentum for the valence band maximum 
and the conduction band minimum must be the same. This type of tunneling is illustrated 
in Eigure 37 (a) and takes place in direct bandgap semiconductors such as GaAs and InP 
as well as indirect bandgap semiconductors where the voltage applied to the junction is 
large enough that the valence band maximum (E point) is in line with the conduction 
band minimum (E point) [7]. 

Conversely, indirect band tunneling occurs when the conduction band minimum 
and the valence band maximum do not occur at the same momentum, as shown in Eigure 
37 (b). As a result, momentum is conserved by scattering agents such as photons or 
impurities and is deemed photon assisted tunneling. Due to the complex nature of photon 
assisted tunneling, it is far less likely to take place when direct tunneling is possible. 
However, it is possible for both direct and indirect tunneling to occur at the same time. 
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Figure 37. (a) Direct band tunneling process, (b) Indirect band tunneling process, 

from [7] 


J. TUNNELING CURRENT 

For quantum tunneling to exist, the electric field in the depletion region must be 
on the order of lO^V/cm. When this is achieved, there is a finite probability that an 
electron will be excited from the valence band directly to the conduction band of the 
semiconductor. 

(U) 

Local band-to-band tunneling models in Silvaco’s ATLAS account for the 
electric field at each node to calculate the generation rate due to tunneling. However, in 
reality the tunneling is non-local and it is therefore necessary to take the spatial profile of 
the energy bands into account as well as the spatial separation of the electrons generated 
in the conduction band from the holes generated in the valence band [9]. To achieve this, 
ATLAS creates a series of 1-D slices through the device to calculate the tunneling 
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current. The sliees are generated as an additional mesh through a series of qty and qtx 
statements, whieh are then plaeed over the top of the original mesh. We will diseuss the 
mesh statements in detail in Seetion V. 

With suffleiently high doping of both the p- and n-region of the junetion, any 
eleetron within the appropriate energy range has a finite probability of tunneling from the 
valenee band to the conduetion band. ATLAS eonsiders eaeh energy level within this 
range and ealeulates the spatial start (-a:,) and end (Xj) point at that level. The tunneling 

eurrent ean then be ealeulated as a funetion of the tunneling probability (T(£■)), the 

effeetive mass (m*), and the potential barrier height (E^ ) as deseribed by Equation 4.2, 

where k is Boltzmann’s eonstant, q the eleetrie eharge, T is temperature, ti is the redueed 
planks eonstant and E is the eleetron energy. 


J(E) 


^ ’ [l + exp[E^^-E]/kr 


EE 


Eqn. 4.2 


ATEAS provides the tools to implement a number of probability equations for 
both loeal and non-loeal tunneling but none proved better at aehieving the forward bias 
than the Wentzel-Kramers-Brillouin method (WKB). 

K. WENTZEL-KRAMERS-BRILLOUIN METHOD (WKB) 

The tunneling probability is given by the WKB method (Eqn. 4.3) in whieh the 
wave veetor (k(x)) of a earner is evaluated as an exponential funetion and integrated 
over the depletion region from -Xj to the Xj. 


Tj □ exp 


-ij " |k(x) dx 


Eqn. 4.3 
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Although there are many ways to implement this equation based on the geometry 

of the potential barrier, the method I will use is based on a triangular geometry (Figure 

( 6 ) 

38) because this is the method used by the ATLAS modeling software. 



Figure 38. Triangular potential barrier for tunneling model, from [7] 


For the triangular potential barrier in Figure 38 the wave vector is given as [7]; 


ki^x) 


2m 


{PE-E) = 


2m 


( 




— - Of'tr 
V2 y 


Eqn. 4.4 


Where PE is the potential energy, E is the incoming electron energy. Eg is the energy 
bandgap, and is the electric field. 


Substituting Equation 4.4 into Equation 4.3 yields [7]; 
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The quantum tunneling probability can be simplified to: 


□ exp- - - Eqn. 4.6 

3 qf^h 


Where is the bandgap of the material, in is the effective mass defined in 
Equation 4.7, q is the electric charge, W is the electric field, and fi is the reduced 
Plank’s constant. 

m* = (xj)(xj) Eqn. 4.7 
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L. SUMMARY 

The tunneling phenomenon presented is based on the probability that an electron 
is able to penetrate a potential barrier if there are occupied and unoccupied states at the 
same energy level on both sides of the potential barrier, the distance between the 
beginning and end of the potential barrier is sufficiently small, and momentum of the 
tunneling electron is conserved. The following Chapter will implement the given 
equations into the Silvaco software. 
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V. SIMULATION SOFTWARE 


A. SILVACO 

Silicon Valley Company (Silvaco) is a leading vendor in technology computer 
aided design (TCAD). Established in 1984 and located in Santa Clara, California, 
Silvaeo has developed a number of exceptional CAD simulation tools to aid in 
semiconduetor proeess and device simulation. Silvaeo has an extensive support team to 
assist with the broad array of semiconduetor technologies. 

B. ATLAS®/DECKBUILD 

The ability to accurately simulate a semiconductor device is eritical to industry 
and research environments. The ATLAS deviee simulator is specifieally designed for 
2D and 3D modeling to inelude eleetrical, optieal and thermal properties within a 
semiconductor device. ATLAS provides an integrated physies-based platform to 
analyze DC, AC and time-domain responses for all semieonductor-based technologies. 
The powerful input syntax allows the user to design any semieonductor device using both 
standard and user-defined material of any size and dimension. ATLAS also offers a 
number of useful deviee examples to assist in your own unique design. 

DeckBuild is the run time environment used to input a eommand file or deck and 

(G) 

is given the extension “.in”. To run ATLAS in the DeckBuild environment, the user 
must first call the ATLAS simulator with the command; 


go atlas 


There is other simulation software available that ean be used in conjunetion with 
DeckBuild, such as ATHENA and DevEdit (Ligure 39), but ATLAS® is the most 
appropriate environment for this researeh. 
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(Run Time Environment) 


Figure 39. Silvaco Inc.’s simulation environment, from [10]. 


C. SYNTAX 

Once ATLAS is called there is a syntax structure that must be followed in order 
for ATLAS to execute the command file successfully. Figure 40 is a list of primary 
group and statement structure specifications. Although there are a few input exceptions, 
the input file statements follow the general format of: 

<STATEMENT> <PARAMETER>=<VALUE> 
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Statements 


Group 


1. structure Specification 


MESH 

REGION 

ELECTRODE 

DOPING 


2. Material Models Specification 


MATERIAL 

MODELS 

CONTACT 

INTERFACE 


3. Numerical Method Selection METHOD 


4. Solution Specification 


LOG 

SOLVE 

LOAD 

SAVE 


5. Results Analysis EXTRACT 

TONYPLOT 


Figure 40. ATLAS command groups with the primary statements in each group, 

from [11]. 

D. STRUCTURE SPECIFICATION 
I. Mesh 

The mesh statement is used to define the structure in an inverted 2D or 3D 
Cartesian grid. The x-axis is positive from left to right, the y-axis is negative from bottom 
to top and the z-axis is positive from the plane of the x-y axis into the workspace. The 
reason for the inverted y-axis is that the manufacturing coordinates are usually described 
as depth below the surface. All coordinates are entered in microns and the spacing is 
used to refine the sharpness and accuracy at an assigned location. ATLAS produces a 
series of triangles to form the mesh based on the user input parameters. Figure 41 shows 
the mesh of a tunnel junction along with the parameters entered in ATLAS . 
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Figure 41. Tunnel junction mesh. 


This particular tunnel junction is only Ipm thick and lOpm wide. The depth (z- 
direction) has a default value of 1pm but was defined as 1 x 10^ to produce an output 
equal to one square centimeter, which corresponds to published data. Building this 
semiconductor as a 1cm x 1cm square would only add nodes to the structure and 
ultimately increase the calculation time. The implementation of the junction to a solar 
cell required a smaller node set which I achieved by limiting the x-y dimensions to 
lOcm^. 

The accuracy of any calculated value in ATLAS is dependent on the refinement 
of the mesh at the junction of interest. The x.mesh is much larger than the y.mesh 
because the variables being measured/calculated are across the junction, which was 
defined horizontally. 
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2. Region 

The region statement is used to separate the initial mesh statement into distinet 
bloeks and sets the initial material parameters that ean be referred to later by region 
number. All meshed areas of a strueture must be assigned a region and the regions must 
be ordered from lowest to highest region. For instance, region 5 cannot be defined before 
region 4. The tu nn el junction is described by the DeckBuild commands along with the 
region lines in Figure 42. 


Data from tjsi.str 



Microns 

Figure 42. Tunnel junction region. 

3. Electrode 

Once the regions are set, the electrodes must be assigned to the desired region so 
that it can be electrically analyzed. The electrode can be assigned to any region or 
portion of a region including the “top” or “bottom” as it is in this research. Because this 
tunnel junction is intended to be used as an interconnect between two solar cells, it was 
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unnecessary to assign the electrodes to metal contacts or large areas that would constitute 
superfluous calculations or error. Therefore, the anode was assigned to the top of the 
junction and cathode was assigned to the bottom. The following statements were used to 
define these parameters: 

telectrodes 

electrode name=anode top 
electrode name=cathode bottom 


4. Doping 

The last required input of the structure specification is the doping statement. The 
doping statement is used to assign the doping level within the previously assigned 
regions. Various properties can be appended to the doping statement to specify how the 
semiconductor was doped and of whether the region is n- or p-type. The desired I-V 
relationship in forward bias will determine the doping characteristics of the 
semiconductor. However, the following statements were developed for a Si junction with 
a forward bias current of approximately 12A at 0.06V: 

doping uniform p.type conc=5el9 y.max=0.5 
doping uniform n.type conc=le20 y.min=0.5 

Here the regions were already defined in a previous statement and were therefore 
referred to by their region number. Region 1 was uniformly doped with a p-type material 
with a concentration of 5el9. Similarly, region 2 was uniformly doped with an n-type 
material at a concentration of le20. 

With the doping statement complete, the structure specifications are complete and 
it is possible to continue with the simulation following the guidelines in Figure 40. 
Unlike the structure specifications, the remaining four groups to be specified do not 
require all statements to run without error. 
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E. QUANTUM TUNNELING MESH 


In order to calculate the quantum tunneling through a potential barrier, finer slices 
of the junction must be taken. ATLAS does this with the definition of the qtx.mesh and 
qty.mesh statements. The direction of tunneling must also be specified with a different 
parameter in the models statement. Figure 43 illustrates the extremely refined qty.mesh 
along with the input parameters to achieve this level of precision. Although the qtx.mesh 
was defined, the slices are only extremely refined in the y-direction because this is the 
direction of tunneling. The new rectangular mesh is superimposed over the original 
triangular mesh and both mesh statements will be evaluated and accounted for in the 
tunneling current. As with all ATLAS spacing parameters, the qty.mesh spacing is 
specified as pm. 


Data from situnnel.str 
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Materials 


Silicon 

Conductor 


qtx.mesh loc=0 spac=l 

qtx. mesh ioc=10 spac=l 

qty. mesh iocs0.3S spac*0.0i 

qty.mesh loc=0.5 spac=0.0005 

qty.mesh loc®0.65 spac=0.0i 


Microns 


Figure 43. Tunnel junction qtx.mesh and qty.mesh. 
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F. MATERIAL MODELS SPECIFICATION 

1. Material 

The material statement relates physical parameters with the materials assigned to 
the mesh. The important material parameters for most standard semiconductors are 
already defined by ATLAS and therefore do not require any changes. We utilized the 
extensive model libraries in ATLAS since the developed tunnel junctions were based in 
Si, Ge, and GaAs which are all very well researched and documented. 

2. Models 

The models statement is essential to the accurate modeling of a particular 
phenomenon because it sets flags for ATLAS to indicate the inclusion of different 
mathematical models, physical mechanisms and other global parameters such as substrate 
temperature. The model statement and parameters used in this research are listed below: 

models temperature=300 srh fermi ni.fermi print qtunn.dir=0 \ 
bbt.nonlocal bbt.forward bbt.nlderivs 


As with all ATLAS statements, a complete list of modeling parameters and their 
definitions can be found in the ATLAS User’s Manual [11]. The most notable 
parameters for the calculation of the quantum tunneling effect are qtunn.dir, bbt.nonlocal, 
bbt.forward and bbt.nlderivs. qtunn.dir specifies the direction of tunneling across the 
junction. Since this device was set up for tunneling across the y-axis, the parameter value 
was set to zero (to specify tunneling in x-direction, qtunn.dir = 1). The bbt.nonlocal, 
bbt.forward and bbt.nlderivs parameters enable the non-local band-to-band tunneling 
model to be used in forward bias and allows the use of the non-local derivatives in the 
Jacobian matrix which helps eliminate convergence problems. 

3. Contact 

Since the tunnel junction is meant to be used as an interconnect of a single crystal 
multi-junction solar cell there was no need to specify any contact parameters. The 
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electrode statements assigning the cathode and anode to the semiconductor were 
sufficient for calculating the current through the junction. 


G. NUMERICAL METHOD SELECTION 

The method statement specifies the numerical methods to be used in the 
calculations specified in the model statement. The quantum tunneling through a junction 
is highly dependent on the carrier concentration. Therefore, it is very important to 
accurately measure the carrier concentration for each iteration. This is done by defining 
the climit parameter in the method statement, which specifies the normalization factor to 
limit the step size of the potential update. The parameter dvmax sets the maximum 
allowed potential updates per iteration. Reducing dvmax dampens the oscillations in the 
small I-V calculations. 

H. SOLUTION SPECIFICATION 

The log/solve/save statements are used to create data files in ATLAS 
simulations. These statements work together to be provide data to be analyzed by other 
functions. This research I have used these three statements as a tool to check bias points, 
energy band diagrams and carrier concentrations as a result of varying stimuli. 

1. Log 

The log statement allows all terminal characteristics generated by a solve 
statement to be saved to a file. 

2. Solve 

The solve statement specifies which bias points are to be applied to produce an 
output. The bias points can be set in a number of different way including step, initial and 
lambda depending on what stimulus is desired. 
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3. 


Save 


The save statement is used to save all node point data into an output file. A 
common use of the log/solve/save statements for the tunnel junction studied is shown 
below: 


log outfile=tjsi.log 

solve name=anode vanode =-l vstep=0.1 vfinal=0.0 
solve name=anode vanode =0.0 vstep=0.01 vfinal=0.2 
solve name=anode vanode =0.2 vstep=0.1 vfinal=1.5 
save outfile=tjsi.log 


Here the output file tjsi.log has been set aside to store data from the solve 
statements and then saved so that the data can be called at a later time. The solve 
statements simply sweep the voltage from -IV to 1.5V in varying voltage steps 
depending on the degree of accuracy at the desired range. With the date stored in an 
outlile, it is ready to be displayed so that it can be analyzed. 

1. RESULTS ANALYSIS 

The tonyplot statement is used to start the graphical post-processor tool. Calling 
tonyplot to visualize tjsi.log displays the output fde in Figure 44. The complete tunnel 
junction input deck can be found in Appendix A. 
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Figure 44. Silicon tunnel junction I-V curve in ATLAS . 
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VI. TUNNEL JUNCTION RESULTS 


After the sueeessful modeling of the basie forward bias tunneling eurrent in Si, it 
was necessary to test the validity of the model. This was accomplished by measuring the 
response of different materials and doping levels to similar bias conditions. 

A. SEMICONDUCTOR MATERIAL 

As an interconnect between two dissimilar solar cells it is important that a tu nn el 
junction can be modeled for materials other than Si. Typical tu nn el junctions used today 
are GaAs and InGaP due to their close lattice matching. The ATLAS input deck 
presented here was developed so that cell parameters are easily changed and the tonyplot 
outputs are easy to compare. 

The doping levels’ initial calculations were based on the Fermi-Dirac equations 
presented in Section II. Figures 45 and 46 illustrate the tunnel junction energy band 
diagrams along with their associated I-V curves for GaAs and InGaP respectively. 



0.1 0.2 0.3 0.4 0.6 0.6 0.7 0.8 0.9 1 
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3 3 

Figure 45. GaAs tunnel junction (p = 3e20/cm n = 9el8/cm) 

{PJcvcv versus volts in a 1pm thick junction). 
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Figure 46. InGaP tunnel junction (p = 2e20/cni^ n = 9el8/cm^) 
{PJcvcv versus volts in a Ipm thick junction). 


By fine-tuning the doping levels, an output current of approximately 200mA was 
realized for both junctions. Though this current is slightly high for a typical current 
produced by a solar cell, we simply solved for 200mA as a baseline. With the successful 
modeling of the two junctions, we were able to manipulate the doping levels and analyze 
the different outputs. 

B. TUNNEL JUNCTION DOPING 

The doping levels of the tunnel junctions were changed for both the n- and p- 
regions to shift the Fermi level and ultimately vary the response to the given bias 
condition. Figures 47 through 52 illustrate the reaction to increased doping of the GaAs 
junction. Figures 53 through 58 show similar results for an InGaP junction. 

Initially, the p-region was varied from a doping level within 3kT (~0.077eV) of 
the valence band, then increased by a reasonable factor (between 2 and 3) to show the 
response in the I-V curve. Pushing the Fermi level further into the valence band 
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corresponds to an increased number of available holes that can recombine with electrons 
tunneling through the junction and therefore increases the current. 

The n-region was varied in a similar manner. However, pushing the Fermi level 
further into the conduction band increases the number of electrons available to tunnel 
through the junction and similarly increases the current. 


1. GaAs Tunnel Junction 


a. p-region variation: 



Figure 47. 


2 

GaAs IV characteristics and energy band diagram (p = le20/cm n = 
9el8/cm^) {AJcvcv versus volts in a 1pm thick junction). 
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Figure 48. GaAs IV characteristics and energy band diagram (p = 3e20/cm^ n 
9el8/cm^) {AJcvcv versus volts in a 1pm thick junction). 



Figure 49. 


GaAs IV characteristics and energy band diagram (p = 6e20/cm n 
9el8/cm^) (A/cm^ versus volts in a 1pm thick junction). 
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b, 


n-region variation 
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Figure 50. 


GaAs IV characteristics and energy band diagram (p = 3e20/cm^ n 
6el8/cm^) {AJcvcv versus volts in a Ipm thick junction). 



Figure 51. GaAs IV characteristics and energy band diagram (p = 3e20/cm n 

9el8/cm^) {AJcvcv versus volts in a 1pm thick junction). 
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Figure 52. 


2 . 


Figure 53. 



3 3 

GaAs IV charaeteristies and energy band diagram (p=3e20/em n=2el9/em ) 
(A/em2 versus volts in a Igm thiek junetion). 


InGaP Tunnel Junction 
a. p-region variation: 
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InGaP IV characteristics and energy band diagram (p = 8eI9/cm n = 

3 2 

9eI8/cm ) (A/cm versus volts in a 1pm thick junction). 
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Figure 54. InGaP IV characteristics and energy band diagram (p = 2e20/cm n 

3 2 

9el8/cm ) (A/cm versus volts in a Ipm thick junction). 



Figure 55. 


■j 

InGaP IV characteristics and energy band diagram (p = 4e20/cm n 
9el8/cm^) {AJcm versus volts in a Ipm thick junction). 
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n-region variation 



Figure 56. InGaP IV characteristics and energy band diagram (p = 2e20/cm n 
6el8/cm^) {AJcvcv versus volts in a 1pm thick junction). 



Data fromtjlnDaP.log 





^ Anode Currer 

(A) 



I ' I I I ' I I I ' I I I ' I I I ' I ' I ' I I I ' I I 

0 0.2 0.4 0.6 0.8 1 1.2 


Distance along line 


Anode Voltage (V!) 


Figure 57. 


InGaP IV characteristics and energy band diagram (p = 2e20/cm n 

3 2 

9el8/cm ) (A/cm versus volts in a 1pm thick junction). 
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Figure 58. InGaP IV characteristics and energy band diagram (p = 2e20/cm n = 
2el9/cm^) {AJcvcv versus volts in a 1pm thick junction). 


C. SUMMARY 

With the characteristic tunnel junction curve successfully modeled, I chose to take 
my research one step further and implement the tunnel junction into a solar cell 
simulation to verify the operation which will be discussed in the next chapter. 
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VII. SOLAR CELL IMPLEMENTATION 


The solar cell implementation was based on research completed by Panayiotis 
Michalopoulos [2] in 2002. Michalopoulos was a student of Prof. Sherif Michael at the 
Naval Postgraduate School and his research was developed using Silvaco software. 

A. SOLAR RADIATION 

The sun’s energy is derived from nuclear fusion reactions where 6 x lO'^kg of 
hydrogen (H 2 ) is converted to helium (He) every second. The mass lost due to the 
conversion is approximately 4 x 10 kg. The energy produced by this loss can be 
determined by Einstein’s famous mass-energy equation (E=mc ) and is found to be 4 x 
10 J. This energy is emitted as electromagnetic radiation primarily at wavelengths from 
the ultraviolet to infrared and radio spectra (0.2 - 3pm) [12]. 

The intensity of the sun’s energy depends on the location of where it is measured. 

2 

Outside the earth’s atmosphere, the intensity is measured to be approximately 1365W/m 
which produces a specific spectral response referred to as air mass zero (AMO). As the 
energy transits the earth’s atmosphere, it is attenuated by ultraviolet absorption due to 
ozone, scattering due to airborne dust, and infrared absorption due to water vapor. AMI 
and AM 1.5 is used for terrestrial applications with incident solar power ranging from 925 
to 963W/ml 

B. I-V CURVE 

As discussed in Chapter II, an electron must release energy to go from an outer 
shell (higher energy state) to an inner shell (lower energy state). Conversely, if a certain 
level of energy is absorbed, such as solar energy, an electron may transit from an inner 
shell to an outer shell or valence shell. In an intrinsic semiconductor this creation of 
carriers leads to an increase in conductivity. However, when a p-n junction is formed in a 
semiconductor the carriers will separate to create electron hole pairs (EHPs). The 
minority carrier will be swept across the junction due to the electrostatic field formed by 
the depletion region and the excess majority carriers created by the photogeneration will 


75 



develop a voltage differential between the two regions. If the junction is then connected 
to a load, current will flow (photocurrent). This phenomenon is known as photovoltaics 
because of the voltage potential generated by the photon (solar) energy. 

As photons generate EHPs in a p-n junction, the standard diode dark current 
characteristic (Figure 59a-blue) is shifted down to produce the illuminated current 
characteristic (Figure 59a-red). Because a solar cell is usually presented as a power 
supply, the illuminated characteristic is most often reflected across the x-axis to show that 
the photocurrent conventionally flows in the opposite direction of the standard diode 
current (Figure 7.1b). 




Figure 59. (a) Dark (blue) and illuminated (red) current characteristic curve of the solar 

cell, (b) Solar cell current characteristic curve with FF. 


The number of FHPs generated in the semiconductor is proportional to the 
photons entering the semiconductor. However, many of the FHPs recombine inside the 
material before it is collected and are therefore not used to produce photocurrent. Thus, 
the maximum obtained power (Pobtamed) than the maximum theoretical power 
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(^theory )• ^hc Tatio of obtamcd power to theoretical power is called fill factor or FF (Eqn. 
7.1). The FF can be seen graphically in Figure 59 (b). 


FF = 


^ obtained 


^ theory 


I V 

max max 

I V 

^sc oc 


Eqn. 7.1 


The power conversion efficiency (n) of a solar cell is given as the ratio of the 
maximum realized power to the incident power (). 


P I V EE 

n = max _ ^SC ^ 

P„„ P„„ 


Eqn. 7.2 


To maximize the power conversion efficiency, the EE, Voc and fie must be 
optimized. 

C. SINGLE CELL 

Initially a single cell for each material in the multi-junction cell was modeled in 
ATLAS . This simulation was based on the code already developed for a Naval 
Postgraduate School thesis by Panayiotis Michalopoulos in March 2002. Both cells were 
illuminated at AMO. The doping levels and thicknesses are illustrated in Eigure 60 and 
62 followed by their corresponding 1-V curve (Eigures 61 and 63). 
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Figure 60. Optimized InGaP cell structure. 
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Figure 61. Optimized InGaP I-V curve. (A/cm versus voltage). 
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Figure 62. Optimized GaAs eell strueture. 


Data from GaAsMigh1.log 



2 

Figure 63. Optimized GaAs I-V eurve. (A/em versus voltage). 
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The InGaP cell achieved a Voc = 1.43V and an Isc = 19niA/cm^. The GaAs cell 

2 

achieved a Voc = I.IV and an Isc = 27.4niA/cm . Both cells were successfully modeled 
with 5% higher power than published data and 2% higher power than the results achieved 
by Michalopoulos [2], 

With the single cells successfully modeled, the tunnel junction was implemented 
into the ATLAS® code. 

D. MULTI CELLS 

Similar to the single solar cells, the multi-junction cell in Figure 64 was 
illuminated at AMO and the I-V curve was recorded (Figure 65). 
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Figure 64. Multi-junction solar cell (InGaP/GaAs) with tunnel junction. 
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Figure 65. Multi-junction solar cell (InGaP/GaAs) I-V curve. (A/cm versus voltage). 


The multi-junction solar cell model achieved a Voc = 2.5V and an Isc = 
18.7mA/cm . As expected, the voltage is the sum of the two cells in the single cell 
models accounting for the small voltage drop across the tu nn el junction. Ideally, both 
cells would generate the same current as in the single cell model and the realized current 
in the multi-junction cell would be the lower of the two. However, because the higher 
energy photons are absorbed by the InGaP top-cell, the GaAs bottom-cell is not able to 
generate as many EHPs and therefore realizes a slightly smaller current than ideal. 

Taking the research one step further, the tunnel junction was implemented in the 
optimized multi-junction cell model. This was done so that the modeled output could be 
compared to empirical data in [3]. The optimized cell has an antireflective coating to 
absorb a greater number of photons and a back surface reflector to reflect some of the 
unutilized photons back into the semiconductor material. These additions help to create 
more EHPs and ultimately a higher current and power density. 

Experimental data for an InGaP/GaAs photovoltaic cell produced a Voc = 2.488V 
and an Isc = 23mA/cm [3]. The optimized multi-junction cell model achieved a Voc = 
2.5V and an Isc = 22mA/cm^ with a EE = 84% and rj = 29% (Eigure 66). This is only 
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3.8% greater power density than the experimental data and demonstrates the ability to 
model a tunnel junetion for use in photovoltaie optimization. 
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Figure 66. 
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VIII. RECOMMENDATIONS AND CONCLUSION 


A. RECOMMENDATIONS 

1. Physical Experiment 

The tunnel junetion model presented here has produeed very promising results. 
For this research the most important characteristic of the tunnel junction is the forward 
bias electron tunneling which was successfully modeled using Silvaco Inc’s ATLAS 
software. Unfortunately, time and resources did not allow for a physical model 
comparison. Although an implemented junction was compared to empirical values, the 
tunnel itself was only compared to theoretical values and responses based on 
mathematical models. 

After further refinement and modeling, it is recommended that a project be funded 
to produce a junction with the precise specifications so that it may be tested and 
compared to the modeled parameters. Developing a more efficient tunnel junction and 
continuing to test the accuracy of the models, will lead to the development of more 
efficient solar cells. 

2. Other Software 

Throughout this study, a number of competing software programs were used in an 
effort to model the forward bias response of a tunnel junction. Some were geared more 
toward other aspects of the semiconductor industries such as MEMs design or computer 
architecture and were therefore difficult to use to implement a tunnel junction into a solar 
cell. Other programs were solar cell centric yet did not have the capacity or tools to 
realize the quantum tunneling effect of the degenerately doped junction. 

There are other software packages available that are comparable to Silvaco’s 
ATLAS and would be beneficial to this research. It is recommended that the utility of 
these software packages be explored and compared. 
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3. Temperature and Radiation Effects 

The use of photovoltaies in spaee applieations is very important. Therefore, the 
tunnel junetion should be adapted and tested under extreme eonditions to ensure the 
aeeuraey of the model. The effeets of temperature and radiation on a tunnel junetion 
should be explored to ensure no anomalies exist. 

B. CONCLUSION 

The value of using simulation software sueh as ATLAS® to model 
semieonduetors has revolutionized the semieonduetor industry. Not only are time and 
money saved by minimizing the need to manufaeture and test ehanges to a design, but 
effleieney is being pushed to the limit. 

It has been have shown throughout this thesis that it is possible to model sueh 
effeets as quantum tunneling in a photovoltaie solar eell. To this point, the tunnel 
junetion in a multi-junetion eell was modeled as a vaeuum to insulate the layers, then 
hardwired to pass the eurrent between eells. Although this is effeetive, the eell being 
modeled is not fully optimized. 

A paper entitled '^Multi-Junction Photovoltaic Model Optimizations for Space and 
Solar Concentrator Applications” [13], was presented by Professor Miehael and the 
author on September 3, 2008, at the 23rd European Photovoltaie Solar Energy 
Conferenee, that was held in Valeneia, Spain. The paper foeused on the applieation of 
multi-junetion photovoltaie eells using tunnel junetion teehnology for use in spaee 
applieations. The tunnel junetion modeling has the potential to push photovoltaie 
teehnology to a new level of optimization and power density and it is neeessary to model 
this phenomenon aeeurately. 
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APPENDIX. ATLAS® SOURCE CODE 


A. MAIN STRUCTURE 

go atlas 

# Definition of constants 

# Mesh 

# X-Mesh 

# Y-Mesh 

# Regions 

# Electrodes 

# Doping 

# QT Mesh 

# Material properties 

# Models 

# Light beams 

# Solving 


B. TUNNEL JUNCTION 

go atlas 

1. X-Mesh 

mesh width=le7 
#x-mesh 

x.mesh loc=0 spac=l 

x. mesh loc=10 spac=l 

2. Y-Mesh 

y. mesh loc=0.0 spac=0.01 

y.mesh loc=0.25 spac=0.001 

y.mesh loc=0.75 spac=0.001 

y.mesh loc=l.0 spac=0.01 

3. Regions 

region num=l material=InGaP x.min=0 x.max=10 y.min=0.5 y.max=l 
region num=2 material=InGaP x.min=0 x.max=10 y.min=0 y.max=0.5 

4. Electrodes 

elec name=anode top 

elec name=cathode bottom 

doping uniform p.type conc=2e20 y.max=0.5 
doping uniform n.type conc=2el9 y.min=0.5 
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5 , 


QT-Mesh 


qtx.mesh loc=0 spac=l 
qtx.mesh loc=10 spac=l 


qty.mesh 
qty.mesh 
qty.mesh 


loc=0.35 
loc=0.5 
loc=0.65 


spac=0.01 
spac=0.0005 
spac=0.01 


6. Models/Output 


models temperature=300 srh fermi ni.fermi print \ 
qtunn.dir=0 bbt.nonlocal bbt.forward bbt.nlderivs 

output band.temp traps u.srh taurn taurp band.param con.band val.band \ 
qfn qfp 

log outfile=tjInGaP.str 
solve init 

save outfile=tjInGaP.str 
tonyplot tjInGaP.str 


method climit=1.0 dvmax=0.05 

solve init 
solve prev 

log outfile=tjInGaP.log 

7. Solve 

solve name=anode vanode =-l vstep=0.1 vfinal=0.0 
solve name=anode vanode =0.0 vstep=0.01 vfinal=0.2 
solve name=anode vanode =0.2 vstep=0.1 vfinal=l.l 
tonyplot tjInGaP.log 

8. Plot 

save outfile=tjInGaP.str 
tonyplot tjInGaP.str 


C. COMMON SECTIONS 
1. X-Mesh 

mesh space.mult=l 

# X-Mesh: surface=500 um2 = 1/200,000 cm2 
x.mesh loc=-250 spac=50 
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x.mesh loc=0 spac=10 
x.mesh loc=250 spac=50 

2. Material Properties 

material TAUN=le-7 TAUP=le-7 COPT=1.5e-10 AUGN=8.3e-32 AUGP=1.8e-31 

# Vacuum 

material material=Vacuum real.index=3.3 imag.index=0 

# Ge 

material material=Ge EG300=0.67 PERMITTIVITY=16 AFFINITY=4 
material material=Ge MUN=3900 MUP=1800 
material material=Ge NC300=1.04el9 NV300=6el8 
material material=Ge index.file=Ge.opt 

# GaAs 

material material=GaAs EG300=1.42 PERMITTIVITY=13.1 AFFINITY=4.07 
material material=GaAs MUN=8800 MUP=400 
material material=GaAs NC300=4.7el7 NV300=7el8 
material material=GaAs index.file=GaAs.opt 

# InGaP 

material material=InGaP EG300=1.9 PERMITTIVITY=11.62 AFFINITY=4.16 
material material=InGaP MUN=1945 MUP=141 
material material=InGaP NC300=1.3e20 NV300=1.28el9 
material material=InGaP index.file=InGaP-l.9.opt 

# AllnP (=InAsP) 

material material=InAsP EG300=2.4 PERMITTIVITY=11.7 AFFINITY=4.2 
material material=InAsP MUN=2291 MUP=142 
material material=InAsP NC300=1.08e20 NV300=1.28el9 
material material=InAsP index.file=AlInP.opt 

# AlInGaP (=InAlAsP) 

material material=InAlAsP EG300=2.4 PERMITTIVITY=11.7 AFFINITY=4.2 
material material=InAlAsP MUN=2150 MUP=141 
material material=InAlAsP NC300=1.2e20 NV300=1.28el9 
material material=InAlAsP index.file=AlInP.opt 

3. Models 

models fermi ni.fermi qtunn.dir=0 bbt.nonlocal bbt.forward bbt.nlderivs 
bbt.std 

4. Light Beams 

beam num=l x.origin=0 y.origin=-5 angle=90 \ 

power.file=AM0silv.spec wavel.start=0.21 wavel.end=4 wavel.num=50 
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D 


INGAAS / GAAS SINGLE AND DUAL JUNCTION CELL 


I. Bottom Cell (GaAs) 

a. Y-Mesh 


# Vacuum 

y.mesh loc=-0.15 spac=0.001 

# Window (0.05 urn) 
y.mesh loc=-0.1 spac=0.01 

# Emitter (0.1 urn) 
y.mesh loc=0 spac=0.01 

# Base (3 urn) 

y.mesh loc=l.5 spac=0.3 
y.mesh loc=3 spac=0.01 

# BSE (0.1 urn) 

y.mesh loc=3.1 spac=0.01 

# Buffer (0.3 urn) 
y.mesh loc=3.4 spac=0.05 

# Substrate (300 urn) 
y.mesh loc=303.4 spac=50 

b. Regions 

# Window AllnP (=InAsP) 

region num=l material=InAsP x.min=-250 x.max=250 y.min=-0.15 y.max=-0.1 

# Emitter 

region num=2 material=GaAs x.min=-250 x.max=250 y.min=-0.1 y.max=0 

# Base 

region num=3 material=GaAs x.min=-250 x.max=250 y.min=0 y.max=3 

# BSE 

region num=4 material=InGaP x.min=-250 x.max=250 y.min=3 y.max=3.1 

# Buffer 

region num=5 material=GaAs x.min=-250 x.max=250 y.min=3.1 y.max=3.4 

# Substrate 

region num=6 material=GaAs x.min=-250 x.max=250 y.min=3.4 y.max=303.4 


c. Electrodes 

electrode name=cathode x.min=-250 x.max=250 y.min=-0.15 y.max=-0.15 
electrode name=anode x.min=-250 x.max=250 y.min=303.4 y.max=303.4 


d. Doping 


# Window 

doping uniform region=l n.type 

# Emitter 

doping uniform region=2 n.type 

# Base 

doping uniform region=3 p.type 

# BSE 

doping uniform region=4 p.type 


conc=lel9 
conc=2el8 
conc=lel7 
conc=2el8 
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# Buffer 

doping uniform region=5 p.type conc=7el8 

# Substrate 

doping uniform region=6 p.type conc=lel9 

e. Solving 


method newton maxtraps=10 itlimit=100 climit=l dvmax=0.05 
solve initial 
solve bl=0.9 

log outf=multiIV.log 

solve vanode=0.0 name=anode vstep=0.025 vfinal=2.7 
extract init infile="multiIV.log" 


extract name="Jsc" 

X.val=0.0 

y. val 

from 

curve(v. 

"anode", 

i."anode") 

where 

extract name="Voc" 
y.val=0.0 

X. val 

from 

curve(v. 

"anode", 

i."anode") 

where 


2. Top Cell (InGaP) 
a. Y-Mesh 


# Vacuum 

y.mesh loc=-0.87 spac=0.003 

# Window (0.03 urn) 

y.mesh loc=-0.84 spac=0.003 

# Emitter (0.05 urn) 
y.mesh loc=-0.79 spac=0.003 

# Base (0.55 urn) 
y.mesh loc=-0.5 spac=0.1 
y.mesh loc=-0.24 spac=0.003 

# BSE (0.03 urn) 

y.mesh loc=-0.21 spac=0.003 

# Buffer (0.03 urn) 

y.mesh loc=-0.18 spac=0.002 

b. Regions 


# Window AllnP (=InAsP) 







region 

0.84 

num=l 

material=InAsP 

X. inin= 

-250 

x.max=250 

y. inin= 

-0.87 

y. max^ 

# Emitter 








region 
0.79 
# Base 

num=2 

material=InGaP 

X. inin= 

-250 

x.max=250 

y. inin= 

-0.84 

y. max^ 

region 
0.24 
# BSE 

num=3 

material=InGaP 

X. inin= 

-250 

x.max=250 

y. inin= 

-0.79 

y. max^ 

region 

0.21 

num=4 

material=InGaP 

X. inin= 

-250 

x.max=250 

y. inin= 

-0.24 

y. max^ 
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# Buffer AllnP (=InAsP) 

region num=5 material=InAsP x.min=-250 x.max=250 y.min=-0.21 y.max=- 
0 . 18 

c. Electrodes 

electrode name=cathode x.min=-250 x.max=250 y.min=-0.87 y.max=-0.87 
electrode name=anode x.min=-250 x.max=250 y.min=-0.18 y.max=-0.18 

d. Doping 

# Window 

doping uniform region=l n.type conc=1.95el8 

# Emitter 

doping uniform region=2 n.type conc=2el8 

# Base 

doping uniform region=3 p.type conc=1.5el7 

# BSE 

doping uniform region=4 p.type conc=2el8 

# Buffer 

doping uniform region=5 p.type conc=0.95el8 

e. Solving 

(1) IscandVoc 

method newton maxtraps=10 itlimit=100 climit=l dvmax=0.05 
solve initial 
solve bl=0.9 

log outf=multiIV.log 

solve vanode=0.0 name=anode vstep=0.025 vfinal=2.7 
extract init infile="multiIV.log" 

extract name="Jsc" y.val from curve(v."anode", i."anode") where 
X.val=0.0 

extract name="Voc" x.val from curve(v."anode", i."anode") where 
y.val=0.0 


3. Stacked Cell 
a. Y-Mesh 

# Vacuum 

y.mesh loc=-0.87 spac=0.003 

# Window (0.03 urn) 

y.mesh loc=-0.84 spac=0.003 

# Emitter (0.05 urn) 
y.mesh loc=-0.79 spac=0.003 

# Base (0.55 urn) 
y.mesh loc=-0.5 spac=0.1 
y.mesh loc=-0.24 spac=0.003 
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# BSF (0.03 urn) 

y.mesh loc=-0.21 spac=0.003 

# Buffer (0.03 urn) 

y.mesh loc=-0.18 spac=0.002 

# Vacuum (0.015 urn) 

y.mesh loc=-0.165 spac=0.002 

# Vacuum (0.015 urn) 
y.mesh loc=-0.15 spac=0.001 

# Window (0.05 urn) 
y.mesh loc=-0.1 spac=0.01 

# Emitter (0.1 urn) 
y.mesh loc=0 spac=0.01 

# Base (3 urn) 

y.mesh loc=l.5 spac=0.3 
y.mesh loc=3 spac=0.01 

# BSF (0.1 urn) 

y.mesh loc=3.1 spac=0.01 

# Buffer (0.3 urn) 
y.mesh loc=3.4 spac=0.05 

# Substrate (300 urn) 
y.mesh loc=303.4 spac=50 

b. Regions 

# Window AllnP (=InAsP) 

region num=l material=lnAsP x.min=-250 x.max=250 y.min=-0.87 y.max=- 
0.84 

# Emitter 

region num=2 material=InGaP x.min=-250 x.max=250 y.min=-0.84 y.max=- 
0.79 

# Base 

region num=3 material=InGaP x.min=-250 x.max=250 y.min=-0.79 y.max=- 
0.24 

# BSF 

region num=4 material=InGaP x.min=-250 x.max=250 y.min=-0.24 y.max=- 

0.21 

# Buffer AllnP (=InAsP) 

region num=5 material=InAsP x.min=-250 x.max=250 y.min=-0.21 y.max=- 
0 . 18 

# Vacuum 

region num=6 material=Vacuum x.min=-250 x.max=250 y.min=-0.18 
y.max=-0.165 

region num=7 material=Vacuum x.min=-250 x.max=250 y.min=-0.165 
y.max=-0.15 

# Window AllnP (=InAsP) 

region num=8 material=InAsP x.min=-250 x.max=250 y.min=-0.15 
y.max=-0.1 

# Emitter 

region num=9 material=GaAs x.min=-250 x.max=250 y.min=-0.1 y.max=0 

# Base 

region num=10 material=GaAs x.min=-250 x.max=250 y.min=0 y.max=3 

# BSF 

region num=ll material=lnGaP x.min=-250 x.max=250 y.min=3 y.max=3.1 

# Buffer 

region num=12 material=GaAs x.min=-250 x.max=250 y.min=3.1 y.max=3.4 
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# Substrate 

region num=13 material=GaAs x.min=-250 x.max=250 y.min=3.4 y.max=303.4 

c. Electrodes 

electrode name=cathode x.min=-250 x.max=250 y.min=-0.87 y.max=-0.87 
electrode name=cc x.min=-250 x.max=250 y.min=-0.18 y.max=-0.18 
electrode name=ee x.min=-250 x.max=250 y.min=-0.15 y.max=-0.15 
electrode name=anode x.min=-250 x.max=250 y.min=303.4 y.max=303.4 

d. Doping 

# Window 

doping uniform region=l n.type conc=1.95el8 

# Emitter 

doping uniform region=2 n.type conc=2el8 

# Base 

doping uniform region=3 p.type conc=1.5el7 

# BSE 

doping uniform region=4 p.type conc=2el8 

# Buffer 

doping uniform region=5 p.type conc=0.95el8 

# Window 

doping uniform region=8 n.type conc=lel9 

# Emitter 

doping uniform region=9 n.type conc=2el8 

# Base 

doping uniform region=10 p.type conc=lel7 

# BSE 

doping uniform region=ll p.type conc=2el8 

# Buffer 

doping uniform region=12 p.type conc=7el8 

# Substrate 

doping uniform region=13 p.type conc=lel9 

e. Solving 

method newton maxtraps=10 itlimit=100 climit=l dvmax=0.05 
solve initial 
solve bl=0.9 

log outf=multiIV.log 

solve vanode=0.0 name=anode vstep=0.025 vfinal=2.7 
extract init infile="multiIV.log" 


extract name="Jsc" 

X.val=0.0 

y. val 

from 

curve(v. 

"anode", 

i."anode") 

where 

extract name="Voc" 
y.val=0.0 

X. val 

from 

curve(v. 

"anode", 

i."anode") 

where 
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4 


Tunnel Junction 


a. Y-Mesh 

y.mesh loc=-0.18 spac=0.002 

# Tunnel emitter (0.015 urn) 
y.mesh loc=-0.165 spac=0.002 

# Tunnel base (0.015 urn) 
y.mesh loc=-0.15 spac=0.001 

b. Regions 

region num=l material=InGaP x.min=-250 x.max=250 y.min=-0.18 y.max= 
0.165 

region num=2 material=InGaP x.min=-250 x.max=250 y.min=-0.165 y.max 
0 . 15 


c. Electrodes 

electrode name=cathode x.min=-250 x.max=250 y.min=-0.18 y.max=-0.18 
electrode name=anode x.min=-250 x.max=250 y.min=-0.15 y.max=-0.15 


d. Doping 

doping uniform region=l p.type conc=8el8 
doping uniform region=2 n.type conc=lel9 


e. QT Mesh 


qtx.mesh 
qtx.mesh 
qtx.mesh 


loc=-250 

loc=0 

loc=250 


spac=50 

spac=10 

spac=50 


qty.mesh loc==-0.165 
qty.mesh loc=-0.1575 
qty.mesh loc=-0.15 


spac=.001 
spac=.0005 
spac=.001 


f. Solving 

method newton maxtraps=10 itlimit=100 climit=l dvmax=0.05 
solve initial 
solve bl=0.9 

log outf=multiIV.log 

solve vanode=0.0 name=anode vstep=0.025 vfinal=2.7 
extract init infile="multiIV.log" 

extract name="Jsc" y.val from curve(v."anode", 1."anode") where 
X.val=0.0 

extract name="Voc" x.val from curve(v."anode", 1."anode") where 
y.val=0.0 
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5 , 


MJ Cell 


a. Y-Mesh 


# Vacuum 

y.mesh loc=-0.87 spac=0.003 

# Window (0.03 urn) 

y.mesh loc=-0.84 spac=0.003 

# Emitter (0.05 urn) 
y.mesh loc=-0.79 spac=0.003 

# Base (0.55 urn) 
y.mesh loc=-0.5 spac=0.1 
y.mesh loc=-0.24 spac=0.003 

# BSE (0.03 urn) 

y.mesh loc=-0.21 spac=0.003 

# Buffer (0.03 urn) 

y.mesh loc=-0.18 spac=0.002 

# Tunnel emitter (0.015 urn) 
y.mesh loc=-0.165 spac=0.002 

# Tunnel base (0.015 urn) 
y.mesh loc=-0.15 spac=0.001 

# Window (0.05 urn) 
y.mesh loc=-0.1 spac=0.01 

# Emitter (0.1 urn) 
y.mesh loc=0 spac=0.01 

# Base (3 urn) 

y.mesh loc=l.5 spac=0.3 
y.mesh loc=3 spac=0.01 

# BSE (0.1 urn) 

y.mesh loc=3.1 spac=0.01 

# Buffer (0.3 urn) 
y.mesh loc=3.4 spac=0.05 

# Substrate (300 urn) 
y.mesh loc=303.4 spac=50 


b. Regions 


# Window AllnP (=InAsP) 

region num=l material=InAsP x.min=-250 x.max=250 
0.84 

# Emitter 

region num=2 material=InGaP x.min=-250 x.max=250 
0.79 

# Base 

region num=3 material=InGaP x.min=-250 x.max=250 
0.24 

# BSE 

region num=4 material=InGaP x.min=-250 x.max=250 

0.21 

# Buffer AllnP (=InAsP) 

region num=5 material=InAsP x.min=-250 x.max=250 
0 . 18 


y.min=-0.87 y.max= 


y.min=-0.84 y.max= 


y.min=-0.79 y.max= 


y.min=-0.24 y.max= 


y.min=-0.21 y.max= 


# Tunnel emitter 
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region num=6 material=InGaP x.min=-250 x.max=250 y.min=-0.18 
y.max=-0.165 

# Tunnel base 

region num=7 material=InGaP x.min=-250 x.max=250 y.min=-0.165 
y.max=-0.15 

# Window AllnP (=InAsP) 

region num=8 material=InAsP x.min=-250 x.max=250 y.min=-0.15 y.max=-0.1 

# Emitter 

region num=9 material=GaAs x.min=-250 x.max=250 y.min=-0.1 y.max=0 

# Base 

region num=10 material=GaAs x.min=-250 x.max=250 y.min=0 y.max=3 

# BSE 

region num=ll material=InGaP x.min=-250 x.max=250 y.min=3 y.max=3.1 

# Buffer 

region num=12 material=GaAs x.min=-250 x.max=250 y.min=3.1 y.max=3.4 

# Substrate 

region num=13 material=GaAs x.min=-250 x.max=250 y.min=3.4 y.max=303.4 

c. Electrodes 

electrode name=cathode x.min=-250 x.max=250 y.min=-0.87 y.max=-0.87 
electrode name=anode x.min=-250 x.max=250 y.min=303.4 y.max=303.4 

d. Doping 

# Window 

doping uniform region=l n.type conc=1.95el8 

# Emitter 

doping uniform region=2 n.type conc=2el8 

# Base 

doping uniform region=3 p.type conc=1.5el7 

# BSE 

doping uniform region=4 p.type conc=2el8 

# Buffer 

doping uniform region=5 p.type conc=0.95el8 

# Tunnel 

doping uniform region=6 p.type conc=8el8 

# Tunnel 

doping uniform region=7 n.type conc=lel9 

# Window 

doping uniform region=8 n.type conc=lel9 

# Emitter 

doping uniform region=9 n.type conc=2el8 

# Base 

doping uniform region=10 p.type conc=lel7 

# BSE 

doping uniform region=ll p.type conc=2el8 

# Buffer 

doping uniform region=12 p.type conc=7el8 

# Substrate 

doping uniform region=13 p.type conc=lel9 

e. Solving 
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method newton maxtraps=10 itlimit=100 climit=l dvmax=0.05 
solve initial 
solve bl=0.9 

log outf=multiIV.log 

solve vanode=0.0 name=anode vstep=0.025 vfinal=2.7 


extract init infile="multiIV.log" 


extract name="Jsc" 

X.val=0.0 

y. val 

from 

curve(v. 

"anode", 

1."anode") 

where 

extract name="Voc" 
y.val=0.0 

X. val 

from 

curve(v. 

"anode", 

1."anode") 

where 
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